Spin-phonon interaction in doped high-Tc superconductors from density functional 

calculations. 
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Effects of coupling between phonon distortions and 
stripe-like spin waves in the CuO plane of HgBa2CuC>4 
are studied by band calculations. Local exchange en- 
hancements depend sensitively on the local structure 
around Cu sites. Interactions where spin waves have 
twice as long wave length as phonon waves can lead to a 
'dip' in the density of states (DOS) below the Fermi en- 
ergy, Ep. This type of interaction is compatible with sev- 
eral seemingly contradictory observations among high-Tc 
cuprates, like the isotope effect, anti- ferromagnetic fluc- 
tuations, stripes and pseudogaps. It can also account for 
a large Tc- 



Numerous studies of high-Tc superconductors have 
revealed unusual properties in the normal and super- 
conducting states. A common feature of the different 
Hg-, Bi-, Y- or La-based high-Tc copper oxides is the 
almost planar CuO sheets, which generates a cylindri- 
cal Fermi-surface (FS) in density functional (DF) band 
calculations. However, real undoped materials become 
anti-ferromagnetic (AF) insulators and it is only within 
a rather narrow range of (usually hole-) doping that 
high-Tc superconductivity can be observed. This del- 
icate balance of the normal state property from being 
metallic or AF insulating leads naturally to the AF pair- 
ing interactional as an alternative to the conventional 
electron-phonon coupling as a possible mechanism for 
high-Tc superconductivity. Complex isotope effects, d- 
wave gap symmetry, magnetic excitations, stripes and 
pseudo gaps all contributeJ.0 the " unconventional" prop- 
erties of these materialsoO. The concept of stripes has 
been developed from models of strong correlation^, where 
Hubbard .model calculations produce gaps in the spectral 
functionstll From DF calculations it has been shown 
that "stripe-like" AF modulations (AFM) of spin- waves 
or phonons ip, the CuO-plane lead to gaps in the one- 
particle DOStl. Such modulations decrease the total en- 
ergy if the gap opens at Ep, which leads to a connection 
between wave length and doping. Taken separately, both 
waves (spin or phonon) give large coupling parameters 
(spin fluctuation or electron-phonon). This motivates a 
closer look at the possible interaction of such modes. We 
show that structural distortions modify local Stoner en- 
hancement and hence the calculations confirm that spin- 



phonon coupling is large. Secondly, it is shown that a 
double gap structure in the DOS occurs when the wave 
lengths of the spin wave is twice that of the phonon mode. 

Interactions between AFM fluctuations and phonon 
modes are modeled for one of the structurally sim- 
plest high-Tc materials, HgBa2Cu04 (HBCO), where 
electronic structure calculations are done for "frozen" 
phonon and spin waves. The DF-calculations use 
the LMTO method based on the local spin density 
approximationllj (LDSA) for supercells of HBCO, in- 
cluding staggered magnetic fields to generate spin-waves 
("stripes") or atomic displacements to model phonons. 
The calculations for the two types of modes are similar, 
but atomic displacements generate multipolar changes of 
the potential and calculations of changes of total energies 
require full potential methods. Long wave length spin 
modulations are simpler since the change in potential is 
spherical within each atom. The calculations are still 
lengthy because of the very large supercells. Therefore, 
in this work we refrain from calculating total energies, 
but we focus on the exchange enhancements and double 
gap structures in thc-|LX)S. The details of the methods 
can be found in refs.M'E3. 

A forced AF (on Cu) configuration in undoped HBCO, 
modeled by LSDA calculations, requires applied field am- 
plitudes of about 23 mRy to open a gap, cf. Table 1. If 
the Ba atoms are moved 0.1 a.u. closer to the CuO plane, 
the critical field is reduced to about 16, while the oppo- 
site movement requires about 30 mRy. Other distortions 
have similar effects, the Ba mode is one example. When 
the atomic distortions are modulated along the cell (but 
with no modulation of the applied AF field) one can no- 
tice that the local AF Cu-moments are largest in a region 
where the Ba are closest to the CuCLplanes. Another 
example is the "half- breathing" modaij, where two pla- 
nar oxygens are moved closer or further away from the 
Cu, so that Cu-sites along [1,0,0] alternatively have two 
O-neighbours at +x and -x closer than normal ("com- 
pressed" Cu), or more distant than normal ("dilated " 
Cu) . The calculated exchange enhancement is roughly 50 
% larger on the dilated than on the compressed Cu. The 
spin induced gap is situated at the stoichiometric band 
filling, 63 valence electrons per HBCO cell or at 1008 
electrons for a 16 formula unit (f.u.) supercell, when the 
AF configuration has no modulation. But as was shown 
in ref.EiJ, if there is a modulation (AFM or phonon) then 
the gap is always found at a band filling corresponding 
to 2 electrons below the position of the gap in the un- 
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modulated case. That is, at 1006 electrons for the 16 f.u. 
cell, at 502 instead of at 504 electrons for the half as long 
cell, etc.. The gain in kinetic energy is largest if the AFM 
wave has the largest moments coinciding with the opti- 
mal local distortions around Cu (Ba close to the CuO 
plane, dilated Cu, etc.). The nodes of the AFM wave, 
where the moments are zero, coincide with the parts of 
the phonon which are most unfavorable to magnetism. 
Furthermore, this implies that the AFM spin wave is 
twice as long as the wave of the phonon modulation. 

As an example we consider a long supercell contain- 
ing 16 f.u., oriented along [1,0,0], cf. fig 1. According 
to the rules for gap positions one expects a gap near 
1006 electrons, caused by one complete spin modulation 
along this cell. We call this gap a 'spin gap'. Simultane- 
ously, there are two wave lengths of a phonon modulation 
within the cell, which gives a gap (called 'phonon gap') 
at a band filling of 1004 electrons. Such features are seen 
in the DOS from the band calculations. The gaps are 
complete with zero DOS in a narrow energy interval, if 
the amplitudes of the atomic distortions or the magnetic 
fields are large. If the amplitudes are moderate, the gaps 
are partial 'pseudo gaps', where the DOS is reduced to 
make 'dips' at the corresponding band fillings. Unpolar- 
ized band calculations for different phonon modulations 
(Ba, apical-O, "breathing-O") show the expected pseudo 
gaps. When the cell contains both the double phonon 
wave and the AFM wave, one should expect two gaps. 
The DOS shown in fig. 2, is calculated for a spin wave 
distribution along the cell which generates a charge mod- 
ulation of two wave lengths along the cell. The strongest 
gap is at Ep (1006 electrons) and it is caused by the spin 
wave. A second smaller 'dip' approximately 2 electrons 
below Ep is caused by a component V q from the charge 
distribution. Note that no phonon distortion is needed 
to see the charge gap in this case, but for a paramagnetic 
calculation with two wave lengths of a phonon, one can 
see a tendency for a gap near 1004 electrons, see fig. 2. 
A similar charge redistribution appears in this separate 
paramagnetic calculation as in the spin polarized case 
with only the AFM wave. There seems to be a mutual 
enforcement of spin and charge gap, as well as increased 
amplitude of the charge transfers in calculations contain- 
ing both the phonon and the AFM wave. 

Another example of this is the half-breathing modes 
within a 8 f.u. cell along [1,0,0]. Two half-breathing 
modes, which are the shortest possible double phonon 
mode that can fit into this cell, induce similar charge 
transfers as the spin wave, as seen in Table 2. The com- 
bination of an unmodulated (AF) spin configuration and 
the phonon wave gives moments that are about 60 per- 
cent larger on the "dilated" than on the "compressed" 
Cu, for the same applied field amplitudes on both types. 
The "compressed" Cu atoms gain charge in the para- 
magnetic calculations, just as the Cu at nonmagnetic 
node positions of a spin wave do. Both the phonon and 
the spin wave enforce the V q coefficient for a (charge) 
gap (at 500 electrons for the 8 f.u. cell). This can in- 



terpreted as a type of valence fluctuation. These have 
been suggested fop,a mechanism for exotic pairing and 
superconductivityt3. The results from the 8 f.u. cell 
show that they are coupled to the spin- phonon modes. 

A few combinations of phonon and spin modulations 
along [1,1,0] show partly similar trends. There are 2 elec- 
trons per supercell separating the gaps induced by the 
spin- and phonon waves. For long wave lengths (long su- 
percells), i.e. for underdoped cases, this means that the 
two gaps should appear relatively close to each other. 
The energy difference between gaps AE w n/N, where 
n and N are the number of states and the DOS per f.u., 
respectively. For instance, for the 16 f.u. cell here, whese 
n is 2/16 per f.u. and N is about 11 states/Ry per f.u.t3, 
the energy difference would be about 140 meV. For over- 
doped cases, when the wave length is short, the two gaps 
should be more separated, although there is some com- 
pensation from the increased DOS as Ep approaches a 
van-Hove singularity at increased-doping. 

Photocmissionllj and tunnelinglJ measurements in the 
Bi-based material reveal a dip at about 80 meV below 
Ep even above Tq- We suggest that this dip corresponds 
to the phonon gap of the DOS, as discussed above. The 
DOS at Ep per CuO plane in the Bi-compound is larscr 
than in HBCO, and increasing further with dopingt3. 
The position of the dip is of the correct order to corre- 
spond to the phonon dip. The origin of this dip was much 
discussed a few years ago, because from tunneling it was 
shown to be pronounced only in the occupied part. This 
is in qualitative agreement with the spin-phonon model. 

Not all phonon displacements produce double gaps. A 
mode with Cu displaced along z makes a relatively small 
difference between the moments on displaced and undis- 
torted atoms. This distortion is symmetric with respect 
to the CuO plane, and the wave length of an AFM mod- 
ulation can be equal to the wave length of the phonon. 
Spin and phonon gaps are weak and appear at the same 
energy. 

As was mentioned above, precise total energies are not 
calculated. However, we can make a number of qualita- 
tive observations of how a strong AF exchange enhance- 
ment and a pseudo gap will change the total energy. The 
spin gap is strongest if a system is close to an AFM tran- 
sition when there is almost no limit for the exchange 
enhancement. Vibrational amplitudes are of the order 
0.1 a.u. at low Tt3. Normal harmonic vibrations follow 
a parabolic potential, where the potential energy U is a 
quadratic function of the displacement x, U — jKx 2 . 
if is a force constant cPU/dx 2 . If, at low T and large 
x, it is possible to gain energy because of a spin wave, 
it will also soften the phonon. The quadratic behavior 
of U will be flattened out (K diminishes with x) which 
allows a further increase of x. Mutual enhancements of 
distortion and moment will enforce both gap structures. 
Several coupled spin-phonon modes may compete at large 
T, when the free energy gain from a gap formation is less 
efficient because of the Fermi-Dirac smearing and ther- 
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mal disorder. A selection at lower T may leave one mode 
with the largest gain in energy to be dominant with a 
more coherent effect on the DOS. 

The Stoner enhancement 5 is useful for a discussion 
of ferro magnetism. An extension for AF can be done 
via generalized spin susceptibility Xq- An approximate 
approach can be based on local 5-factors (defined within 
each site) calculated for AFM waves. Here they are de- 
fined as the exchange splitting of the logarithmic deriva- 
tive of the Cu-d band divided with the applied field (e) at 
the same site. The local Stoner factor 5 defines the local 
exchange integral I, S = NI = 1 — 1/5. An important 
decrease of the total energy is coming from the opening 
of a gap via the potential component V q , which in the 
idealized case is close to — ^N q V 2 , where V q is the band 
shift at q and N q is a fraction of the total DOS that is af- 
fected by the splitting. The total energy has also kinetic 
and exchange contributions from each site, which can be 
summed up, if each site is independent of its neighbors. 

E = £(A l£ 2 - N 2 I t e 2 ) = (1) 

where m = Ne is the magnetic moment, and it is assumed 
that N and 5 are not much varying from site to site. 
The additional hybridization energy due to nonequivalent 
neighboring atoms in the AF case is accounted for in local 
5-values, defined from spin-polarized results for the real 
spin-wave configuration. All these contributions make 
5 to vary non-linearly with the field, as can be seen in 
Table 1 for the simple AF configuration of the undoped 2 
f.u. cell. It is noted that the LSDA potential is too weak 
to explain a spontaneous AF order and a full gap. The 
local / (or N) needs to be almost twice as large to get 
a gap. A similar reasoning can be made for the doped 
cases, with modulations of AFM waves. 

Bardeen-Cooper-Schrieffer (BCS) theory without 
strong coupling effects is often sufficient to understand 
the mechanism for superconductpribz. based on electron- 
phonon or paramagnon couplingEjo. A potential per- 
turbation, as from a phonon, with Fourier component 
V q leads to a "Cooper" pairing of electrons of opposite 
momentum ±q. The BCS relation for Tc contains the 
electron-phonon coupling A p h; 

k B T c = 1.13«lueap(-l/Aph) (2) 

where Qd is a characteristic phonon frequency. Here the 
coupling is for electronic states of opposite spin. Since 
Xph ~ NV 2 /K = NV 2 /MQ 2 D , where V is matrix el- 
ement of the perturbation potential due to a displace- 
ment, there is normally low coupling for stiff phonons. 
High phonon energies in combination with fairly large 
Aph, estimated from band theory, can hardly ex plain the 
very large observed Tc's in the copper oxidestall 2 ], al- 
though the,cO|le of phonons on superconductivity cannot 
be ignoxedHiiil. An analogous formulation for equal spin 
pairingEa is; 

k B Tc = 1.13fia; s /exp(-l/A s f) (3) 



The prefactor u> s f is the frequency of the spin fluctuation 
and A s f is the coupling strength. Tii^e properties are in 
the ferro magnetic case related astSEil: 

hw sf = 2/NS ; A sf = \SS 2 (4) 

As for phonons there is an antagonistic dependence of 
Tc on frequency and coupling strength. For instance, if 
the spin coupling is large via a large 5-factor there is a 
softening, so that u> s f and T c goes down. The equal spin 
pairing has been suggested to gi\p,small T c 's near the 
FM transition in C15 compoundsti Extreme softening 
leads to static spin for magnetic coupling or to a struc- 
tural transformation for phonon coupling. The major 
contribution to the matrix element V for an AFM per- 
turbation comes from the opening of a gap. An approx- 
imate value can be obtained from the gap model, where 
the gap is equal to the maximum value of the exchange 
enhanced field. The model permits for an evaluation of 
the matrix element V ~ 1, which with the denominator 
K = d 2 E/de 2 = 2N/S makes A sf « ±5. Eq. | is for- 
mally as for ferro magnetic cases, whereby one can use 
huj s f ~ 2/NS. Near the region for AF fluctuations, one 
can expect very large A s f , but unless a there is a mecha- 
nism for a large w s f, it will not produce a large Tc- 

As was discussed above, Cu-moments can vary 50 per- 
cent or more along a phonon with realistic distortions. 
This implies strong nonlinear variations as function of 
time and position, between low and high values of local 
Stoner enhancements and uJ s f- If the development of a 
spin wave is strongly coupled to the atomic distortions, as 
is suggested from the calculations, it is unlikely that the 
excitation spectrum for the spins would be given by u) s j , 
since that quantity is independent of the phonons. It is 
more realistic if the spin excitation spectrum is connected 
to The nonlinear variations of 5 permits a rapid de- 
velopment of a spin wave (ui s f — * &>d) when the local 5 
is small, alternated by strong coupling when 5 is large. 
Simultaneous large of values of the prefactor and 5 can 
explain a much larger Tc's than if Tc is evaluated from 
spin fluctuations alone. A major difficulty for quantita- 
tive results from LSDA calculations is that the tendency 
towards AFM is too weak. As indicated above, the ex- 
change integral (or N) needs to be larger to explain a gap 
in the undoped system. If such an increase can be justi- 
fied for the modulated cases, it will make the difference 
between high and low 5 values along a phonon distortion 
even more spectacular. The average LSDA values for 5 
are about 0.5 in the modulated cases, corresponding to a 
A s f of the order 1. Omission of effects from strong cou- 
pling and non-linearity in a BCS-like formula with a fre- 
quency prefactor of 50 meV as for oxygen phonon modes, 
leads to a large Tc of the order 200 K. This is very ap- 
proximate, but it shows qualitatively the possibilities of 
having a large Tc from large AFM enhancements in com- 
bination with high excitation frequencies. Impurities and 
strong electron phonon coupling are destructive to su- 
perconductivity based on the paramagnon mechanismlia. 
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However, with spin-phonon coupling there are mutual en- 
forcements between phonons and spin waves which can 
make Tq more insensitive to perturbations or disorder. 
Spin-phonon coupling would be consistent with a num- 
ber of unusual observations among the high-Tc materials 
such as an isotope effect on Tc and pseudo gaps, even 
in the case of pairing mediated by spin fluctuations. A 
crucial test of spin-phonon coupling would be to measure 
the isotope effect on the spin excitation spectrum, since 
large isotope shifts are expected from the model. 

The calculations show that the conditions for AF de- 
pend sensitively on the local structure around Cu. This 
effect is probably more important than in the present 
results, since LSDA calculations underestimate the ten- 
dency towards AF. Therefore, it can be concluded that 
spin-phonon coupling is important in these materials. 
Similar charge transfers caused by phonons and spin 
waves, or charge transfers induced by stripe-like spin 
waves alone, give a second gap in the DOS just below 
the main gap at Ep. The energy difference between the 
gaps depends on doping, and as was shown previously^, 
the wave length of the stripes are predicted to vary with 
doping as well. Other effects from spin-phonon coupling 
are consistent with many observations in the high- Tc ox- 
ides. In particular, this coupling might be important for 
the mechanism of high-Tc superconductivity. 
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TABLE I. Calculated local Stoner enhancement S on Cu, 
gap in the DOS E g (mRy), as function of the applied field 
(mRy) on Cu in undoped AFM Hg2Ba4Cu208 for normal 
structure and a distorted structure where Ba is moved 0.1 
a.u. towards the CuO plane. 
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TABLE II. Charges and moments on the two types of Cu 
in a 8 f.u. "half-breathing" modulation along [1,0,0]. The 
O-displacement in the phonon is 0.1 a.u.. Cu(l) are "dilated" 
and Cu(2) "compressed" sites, see the text. "AF" has no 
modulation of the anti- ferromagnetic configuration. "AFM" 
is with modulation. 
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Moment 
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10.32 
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10.40 




AF+phon. 
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0.34 
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AFM+phon. 
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FIG. 1. The CuO plane for a supercell of 16 f.u. along 
[1,0,0] with modulated AFM moments. The filled circles are 
Cu and open circles the plane oxygens. The arrows indicate 
up- and down-moments on Cu. 
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FIG. 2. The total DOS for a 16 f.u. cell in fig. 1. The thin 
line is normal paramagnetic, the heavy line is with an AFM 
modulation, which makes a perturbation of the Fourier coef- 
ficient V q in the spin potential. Ef is at 1006 electrons/cell 
(0.125 holes per f.u.). A second "dip" in the DOS is found 
at about 2 electrons/cell below Ef- The origin of this dip 
is found in a modulation of the charge (unpolarized Cu have 
more charge than Cu with moments) which gives a potential 
perturbation in the coefficient Vq, where Q = 2 • q. The bro- 
ken line is for a paramagnetic calculation, where two wave 
lengths of "half- breathing" distortions are introduced along 
the cell, see the text and table 2. 
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